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ABSTRACT 



Aims. We study the blazar nature of the high-redshift Flat-Spectrum Radio Quasar PKS 2149 - 306 (z = 2.345) by investigating its long-term 
behavior. 

Methods. We analyzed all publicly available optical-to-X-ray observations performed by XMM-Newton, Swift, and INTEGRAL. 
Results. PKS 2149 - 306 is one of four blazars at z > 2 that have been observed in the hard-X-ray regime with both the BAT and ISGRI 
instruments. Observations acquired almost 1 year apart in the 60 - 300 keV energy band in the object rest frame, exhibit no noticeable change 
in spectral slope associated with a flux variation of more than a factor of two. Swift data appear to show a roll-off below ~ 1 keV, which 
becomes increasingly evident during a ~ 3-day time-frame, that can be explained as the natural spectral break caused by the Inverse Compton 
onset. The broad-band spectra allow us to identify two different states. The SED modeling suggests that they can be interpreted by only a 
change in the bulk Lorentz factor of the jet. 
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1. Introduction 

PKS 2149 - 306 is a Flat Spectrum Radio Quasar (FSRQ) 
located at z — 2.345 (Wilkes 1986). In the radio maps, it 
exhibits a compact but non-point-like structure (e.g., Ojha 
et al. 2005) of « 1 Jy flux density over a wide radio fre- 
quency range (0.8 - 8.4 GHzfl The few observations of the 
optical counterpart show no significant changes in all filters 
with my ~ 18 (e.g., Francis et al. 2000). The object hosts a 
bright X-ra y source extensively observed in the p ast by ROSAT 
and ASCA dSiebert et al. 19961: ICappi et al. 19971) . BeppoSAX 
( Elvis et al. 2000l). Chandra dFang et al. 200 lb . XMM-Newton 
(Ferrero & Brinkmann 2003), and more recently by Swift 
dSambruna et al. 20071) . In the X-ray domain, the source ex- 
hibits a hard spectrum with variable photon index and flux. 
From early ASCA observations, the presence of an emission 
line at * 17 keV in the blazar frame was interpreted as a highly- 
blues hifted Fe Kg fine produced by an outflow at speed v ~ 
0.7c dYaqoob et al. 19991) : however, Chandra high-resolution 



1 Data from NASA/IPAC Extragalactic Database 
http : / /nedwww . ipac . caltech . edu/ index . html 



(NED) 



spectroscopy did not confirm this finding dFang et al. 2001b . 
Page et al. (2004) found no evidence for a disk/torus Fe Ka 
emission line. 

As for other high-z FSRQs, a controversial low-energy 
photon deficit has bee n claimed by some authors but not 



confirmed by oth ers ( Siebert et al. 1996t ICappi et al. 1997 



Elvis et al. 2000: 



l Yaqoobetal. 19991: ISambruna et al. 20071 
iFerrero & Brin kmann 2 0031: IPage et al. 2005b . This putative 
deficit was usually ascribed to an absorbing cloud in the quasar 
reference frame and modeled by additional intrinsic absorption 
in addition to the Galactic one. However, considering the blazar 
nature of PKS 2149 - 306, it is difficult to explain how an ab- 
sorbing cloud can survive in front of a jet. On the other hand, 
interpreting this flux deficit as a natural break in the Spectral 
Energy Distribution (SED) due to the low-energy tail of the 
electron population leads to a coherent explanation of the ob- 
served spectral changes in PKS 2149 - 306 (cf. the case of 
RBS 315, Tavecchio et al. 2007). 

PKS 2149 - 306 is one of the few FSR Qs observed in th e 
hard-X-ray regime after the BeppoSAX era (Elvis et al. 2000b . 
A recent Swift observation dSambruna et al. 20071) confirmed 
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that the source displays a Compton-dominated spectrum, typ- 
ical of FSRQs, with a Compton peak probably located at ap- 
proximately hundreds of keV (observer frame). 

In this work, we present a serendipitous detection of 
PKS 2149-306 in INTEGRAL archival data To date, the source 
catalog of the INTEGRAL/lSGRl detector dLebrun et al. 2003b 
includes 18 blazars, among which 6 are located at z > 2 
(e.g., Bird et al. 2007). The AGN catalog of the Swift/BAT 
instrument dBarthelmv et al. 2005b contains comparable num - 
bers, with 6 high-redshift objects dBaumgartner et al~20 08). 
Including PKS 2149 - 306, there are only four objects detected 
by both ISGRI and BAT. We performed a spectral analysis of 
INTEGRAL data reanalyzed all publicly available data for the 
source from the XMM-Newton and Swift archives. Here, we dis- 
cuss the broad-band spectral variation of the source over sev- 
eral years. 

2. Data analysis 

Table Q] lists all observations of PKS 2149 - 306 analyzed in 
this paper, including starting dates and exposure times. For all 
observatories, data were processed using the latest versions of 
the specific software and calibration files. The output spectra 
of XMM-Newton/EPIC and SwiftPiRT were rebinned to con- 
tain at least 30 counts per bin. Spectral analysis was carried out 
by means of xspec v. 12 . 4. 0. In all X-ray model fits, we in- 
cluded Galactic extinctio n and fixed the colum n density to be 



N a = 1.63 x 10 2U cm" 2 dKalberla et al. 2005h . Our fit results 
are given in Table [2] 



2.1. INTEGRAL 

We achieved a serendipitous detection of PKS 2149 - 306 us- 
ing t he INTEGRAL public archive data; the source was in the 
IBIS dUbertini et al. 20031) field of view during November 2004 
(Revolutions 254, 258 - 260), while pointing to NGC 7172. 
Data were analyzed with the Off-line Sc ientific analysis 
Software OSA 7 . dCourvo i sier et al. 2003 ). The imaging 
analysis of IBIS/ISGRI data dGoldwurm et al. 20031) provided 
a source detection of signal-to-noise ratio S/N= 4 in the 
20 - 100 keV energy band with exposure of 366 ks. The 
faint detections in single pointings did not allow us to inspect 
light curve variations. The total spectrum was extracted using 
a rebinned response matrix with 5 energy bins in the energy 
range 20 - 140 keV. The ISGRI spectrum can be modeled by 
a single power-law with photon index F = 1.5*2? and flux 
(2.1 + 0.2) x lO" 11 erg cnr 2 s _1 in the 20 - 100 keV energy 
band (Table©. 

The source is below the capabilities of both JEM-X and 
PICsIT instruments, with upper limits of 10~" erg cirT 2 s _1 in 
8-14 keV and 3.5 x 10~ 10 erg cirT 2 s _1 in 252 - 336 keV 
respectively. 

SPI data were not considered since the source NGC 7172, 
the primary target of the observations, falls on the border of the 
instrument Point-Spread Function (PSF) (2.5 degrees) centered 
on PKS 2149-306. 

OMC data were not included in the analysis because a con- 
taminating source is located within the instrumental 25 arcsec 



PSF FWHM of PKS 2 149 - 306 (see the SwiftfUVOT data anal- 
ysis below). 

2.2. XMM-Newton 

XMM-Newton observed PKS 2149 - 306 on May 1, 
2001, starting at 10:53 UT (ObsID 0103060401). We 
reanalyzed these data, already presen ted by Ferrero & 



Brinkmann ( 2003). Data fro m EPIC-PN (Struderet al. 2001), 
EPIC-MOS dTurner et al. 200 lb and OM dMason et al. 200 lb 
were processed, screened, and analyzed using the same pro- 
cedure described in Foschini et al. (2006), but with the SAS v 
7.1.0 software and calibration file release of July 16, 2007. 
No periods of high-background were detected, and the net ex- 
posures with the individual detectors were 20.3 ks with PN and 
24 ks with MOS1 and MOS2. 

The EPIC-PN light curve (not shown here) displays no sig- 
nificant flux changes (x 1 probability of constancy 0.43) with an 
upper limit to the fractional variability < 3% (3cr). 

Our spectral analysis confirms the results of Ferrero & 
Brinkmann (2003) that the best-fit model was a single power- 
law function with parameters reported in Table [2] The normal- 
ization constants between the MOS instruments and PN are 
consistent with unity (Cmosi-pn = 1.00 ± 0.02, Cmosi-pn = 
1.03 ± 0.02). In panel a of Fig. Q] the residuals of the power- 
law model are plotted in units of <x and with error bars of 
size 1 cr; for a clearer visualization, only PN data are shown. 
We applied more complex models: a power-law function with 
Galactic extinction and an absorption excess in the source 
frame (wabs*zwabs*zpo) and an absorbed broken power-law 
function (wabs*bknpo). Both provided comparable values of 
X 2 r , but the parameters were not well-constrained or they con- 
verged to unphysical quantities. 

We also checked for the presence of neutral iron line emis- 
sion in the quasar frame but we found no indication of such a 
feature (confirming Page et al. 2004). 

The OM observed the blazar for ~ 1000 s with UVW2 
filter (centered at 212 nm): no detection was found, with a 
magnitude lower limit of 19 (3cr), corresponding to a flux 
< 3.4 x 10~ 13 erg cirT 2 s" 1 . 



2.3. Swift 

Swift observed PK S 2149 - 306 twice in 2005 
dSambruna et al. 2007b : ObsID 00035242001 started on 
December 10, 2005 00:47 UT, and ObsID 00035242002 
started on December 13, 2005 20:13 UT with exposure 
times of 3.1 ks and 2.3 ks, respectively. For the reduction 



and analysis of the data from the three instruments onboard 
the Swift satellite, we used the HEASoft v. 6.4 package, 
together with the CALDB updated on Feb ruary 13, 2008. 

The data from the BAT instrument dBarthelmv et al. 2 005), 
optimized for the 15 - 195 keV energy band, were binned, 
cleaned of hot pixels, and deconvolved. The intensity maps 
from the individual pointings were integrated by using the 
corresponding variance maps as a weighting factor, and the 
resulting exposure simultaneous to the XRT and UVOT ob- 



Table 1. Observation Log. 
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Instrument 


ObsID 


Start Date 


Exposure 






[YYYY-MM-DD] 


[ks] 


INTEGRA L/ISGRI 


Rev. 254, 258 - 260 


2004 - 1 1 - 1 1 


366 


XMM/EPIC 


0103060401 


2001 - 05 - 01 


20 


SwiftfKKF (OBS. 1) 


00035242001 


2005 - 12 - 10 


3.1 


Swift/XKT (OBS. 2) 


00035242002 


2005 - 12 - 13 


2.3 


SwiftfBAT 


9-month Survey 


mid Dec 2005 


2800 



For each instrument, the table indicates the observation ID or the revolution number, the starting date of the observation, and the effective 
exposure for PKS 2149 - 306. BAT data refer to the 9-month survey, starting in mid-December 2005. BAT exposure time is computed by 
integrating over all pointings with the source PKS 2149 - 306 within the FOV of the BAT instrument. 



servations was 5.4 ks. No source was detected; the corre- 



sponding 3<x upper limits were 2.7 x 10 10 erg cm 2 s 1 and 



,-2 „-l 



4.3 x 1(T 1U erg cirT 2 s _l in the 20-40 keV and 40-100 keV en- 
ergy band, respectivel y. However, the bla zar was detected in the 
9-month data survey (ITueller et al. 2008), and the source spec - 
trum is publicly available online (rBaumgartner et al. 2008 ). 
The fitting of a power-law model is optimal for a photon in- 
dex of T = 1.5 + 0.4 and a flux of 4.8 x 10~ n erg cirT 2 s -1 in 
the 20 - 100 keV energy band. 

The X-Ray Telescope (XRT, Burrows et al. 2005), operat- 
ing in the 0.2 - 10 keV energy band, was set to work in photon 
counting mode during the two pointings on PKS 2149 - 306. 
Data were processed and screened by using the xrtpipeline 
task with the single to quadruple pixels events selected (grades 
0-12). 

The spectral analysis of the XRT data was presented in 
Sambruna et al. (2007) as a joint fit of the two observations. 
Here, we consider the two data sets separately and we fit each 
with three models, all including Galactic extinction: a power- 
law (wabs *zpo), a power-law with an extra-absorption compo- 
nent in the source rest frame (wabs*zwabs*zpo) and a broken 
power-law (wabs*bknpo). 

The spectrum of the first observation (2005-12-10) is well 
described by a single power-law with Galactic absorption 
(Table|2]i. Residuals are shown in panel b of Fig.[T](in units of 
1 cr). Although the broken power-law and the power-law with 
extra absorption models yield ^ = 1 .2 (for 27 and 28 dof, re- 
spectively), the spectral parameters are not well constrained or 
they converge to unphysical values. For the second observation 
(2005-12-13), we found that the broken power-law function 
provided the best-fit model, although both the power-law model 
and the power-law function with an extra-absorption compo- 
nent described the observational data accurately. In Table|2]we 
present the best-fit parameters for the three applied models. The 
residual plot of the power-law model (Fig.Q] panel c) shows a 
hint of curvature, which is not present when the broken power- 
law model is assumed (Fig. Q] panel d). The curved model is 
also supported by the F-test analysis, showing that the flatter 
power-law function below ~ 1.7 ke V is required with a prob- 
ability of 98.6%. The quality of data acquired to date does not 
allow us to distinguish firmly between the intrinsically curva- 
ture (broken power-law) and the extra-absorption model. 



UVOT dRoming et al. 20051) data analysis was performed 
by the uvotmaghist task using a source region of 5" radius 



for the optical and 10" for the UV filters. Since a contami- 
nating source is located close to PKS 2149 - 306, the back- 
ground was evaluated from a nearby region of 60" radius. 
The observed magnitudes were computed as the mean values 
of recorded data, since they do not exhibit significant varia- 
tions: V 54 3nm = 17.4 ± 0.2, B 4 34nm = 17.7 + 0.1, f/ 3 44nm = 

17.1 + 0.1, £/VWl 2 9inm = 17.9 + 0.2, <yVM2 2 3inm > 19.1, and 
[/VW22i2nm ^ 19.4. The optical magnitudes are in accordance 
with values observed by Francis et al. (2000). With respect to 
the previous analysis of Sambruna et al. (2007), we obtain com- 
parable magnitudes for all filters but UVM2 and UVW2, for 
which we found an upper limit. The discrepancy is probably 
due to the use of updated calibration files. 

2.4. Joint fits 

We attempted a joint fit of the spectra provided by EPIC and 
ISGRI data. This is justified by the fact that the source shows a 
low flux level in both ISGRI and EPIC data and that the ISGRI 
data fall on the extrapolation of the EPIC spectrum (the normal- 
ization of the ISGRI spectrum with respect to the EPIC spec- 
trum is C = 1.0 ± 0.4). We then argue that the blazar is in the 
same spectral state during the two observations, despite being 
separated by ~ 4 years. The joint fitting procedure found a sin- 
gle power-law function to be the best-fit model with parameters 
given in Table [2] 

We performed a joint fit to the XRT observations and the 
BAT spectrum. Best-fit models are given in Table [2] We note 
that for the observation of 2005-12-13 the normalization of 
BAT with respect to XRT is higher than unity (C = 2.3+^): 
this probably reflects the flux variability of the source, given 
the fact that the BAT spectrum relates to data integrated over 
a period of 9 months. Since the XRT spectrum was retained 
as reference, we note that the model flux is comparable to that 
obtained for the EPIC+ISGRI spectrum. 

3. Discussion 

3.1. Spectral variation 

The present analysis, and past results, shows that the hard-X- 
ray emission in the ~ 60 - 300 keV interval (object rest frame) 
of PKS 2149 - 306 does not vary in spectral shape: indeed 
spectra obtained almost a year apart by INTEGRALflSGRl and 
Swift/BNI show a constant (within errors) photon index, which 
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Table 2. Fit results. 
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X-rays (0.2 - 


10 keV) 






Instrument 


Model 


T or F\ 


Additional Parameters 


^ 2-10 keV 





XMMfEPIC wabs*zpo 1.42 ±0.01 0.56 ± 0.01 1.08/977 

Swift/XKT (OBS. 1) wabs*zpo 1.44 ±0.07 1.2 ±0.1 1.24/29 

Swift/XKT (OBS. 2) wabs*zpo 1.38 ± 0.09 1-04^m 1.63/15 

wabs*zwabs*zpo 1-55+° jj? N. = 10 22 cnr 2 °- 97 -oo8 1-42/14 

wabs*bknpo 0.98^ 2 « £„reak = 1.9^ keV r 2 = 1.86^ 0.8 ±0.1 0.98/13 



Hard X-rays (20 - 100 keV) 

Instrument Model T ^20-100 kcv ;tr/dof 

INTEGRALflSGRI zpo 1.5^ 2.1 ±0.2 0.51/3 

SwiftfB.NI zpo 1.5 ±0.4 4.8 ± 0.5 0.48/2 



Joint fit (0.2 - 100 keV) 

Instrument Model r or Ti Additional Parameters ^2-100 kcv ^ 2 /dof 

EPIC+ISGRI wabs*zpo 1.42 ±0.01 1.08/980 

XRT(OBS. 1) + BAT wabs*zpo 1.44 ± 0.07 6.4+g-| 1.2/32 

XRT (OBS. 2) + BAT wabs*bknpo 0.9^ £ break = 1.6+{^ keV f 2 = 1.7+°| 3±1 0.92/16 



Models and parameters for single instruments and joint fit. Applied models are: a power-law in the source rest frame (zpo), a power-law with an 
absorption component in the object rest fram e (zwabs*zpo) and a broken power-law (bknpo). All models include Galactic absorption (wabs) 
with column density N H = 1.63 x 10 20 cirr 2 dKalberla et al. 20051) . Flux is given in units of 10 -11 erg cirr 2 s _I . 



is also consistent with the slope previously inferred for data 
acquired by BeppoSAX in 1997 (r = 1.37 + 0.04 - Elvis et 
al. 2000). We recall that the spectra refer to an average state 
covering ~ 10 days for ISGRI and ~ 9 months for BAT. 
Despite the lack of spectral variability, the integrated flux in 
the 20 - 100 keV band exhibits a variation of more than a fac- 
tor of two from the lowest value of ISGRI to the highest value 
for BAT. A remarkable change in the hard-X-ray flux is not 
uncommon for blazars, and it was observed, for instance, in 
Swift J1656.3 - 330 2, although with a different spectral slope 
dMasetti et al. 20071) . 



hand, in the blazar scenario the soft photon deficit and its 
variation can be naturally interpreted a s a spectral break due 



In the 0.2 - 10 keV energy band, our data analysis sug- 
gests a more complex behavior. As discussed in Ferrero & 
Brinkmann (2003), the source spectrum acquired in 2001 
(XMM-Newton/EPIC data) is well described by a single power- 
law function with little evidence that a more complex model 
is necessary. Instead, the two observations with XRT in 2005 
suggest a spectral change from a power-law model in the first 
observation (2005-12-10), to a curved model in the second 
pointing (2005-12-13), with a hint of curvature developing over 
~ 3 days. As for other high-z radio-loud quasars, this feature 
is usually interpreted as an extra-absorption component. This 
absorbing cloud should be associated with the quasar frame, 
since no Damped Lya features reveal absorption by interven- 
ing material. Although short-term variability due to changes 
in the intrinsic abs orption was observed in Seyfert galaxies 
(Ris aliti et al. 20071) . a variation on time intervals of days is 
difficult to reconcile with the blazar picture since the absorb- 
ing cloud should be located at a distance of ~ 10 15 cm from 
the black hole, below the length scale on which jet dissipa- 
tion is supposed to occur in this class of objects. On the other 



to the Inverse Comp ton regime onset ( Tavecchio et al. 2007 



Sambruna et al. 2007b . Indeed a low-energy roll-off appears 
to be a common property of FSRQs, for distant objects such 
as MG3 J225155 + 2217 (z = 3.668, Maraschi et al. 2008, 
RBS 315 (z = 2.69, Tavecchio et al. 2007), and nearby FSRQs 
such as 3C 454.3 (z = 0.859, Ghisellini et al. 2007). For high-z 
FSRQs, the break signature falls in the soft X-ray band, where 
this feature can be blurred by contribution from the primary 
seed photon inp ut and the bulk C omptonization of cold elec- 
trons in the jet dCelotti et al. 20071) . Apart from spectral anal- 
ysis, a variability study of simultaneous data in the optical to 
X-ray bands can provide clues to the or igin of the intrinsic cur- 
vature in the continuum (IFoschini 2 008). 



3.2. SED modeling 

We define two different states of PKS 2149 - 306: a "low- 
flux" state given by XMM-Newton and INTEGRAL points and 
a "high-flux" state studied with XRT and BAT data, for which 
we obtained a high flux level separately in the [2 - 10] and 
[20 - 100] keV band. We chose the second observation of 
XRT (on 2005-12-13), since the flux is comparable to the first 
one but exhibits spectral bending at low energies. In Figj2] 
the SEDs for both states are depicted separately in blue and 
red. Optical/UV data were corrected for Galactic absorption 
according to Cardelli et al. (1989) and using the zero points 
given in the calibration files. X-ray points are corrected only 
for the Galactic absorption. The multiwavelength spectra are 
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Energy (keV) 



Energy (keV) 



2 5 
Energy (keV) 



Energy (keV) 

Fig. 1. Residual plots in units of cr with error bars of size 1 
sigma. Panel a: XMM-Newton data with respect to a power- 
law model in the object rest frame with Galactic absorp- 
tion; Panel b: XRT (2005-12-10) with a single power-law; 
Panel c: XRT (2005-12-13) with a single power-law; Panel d: 
XRT (2005-12-13) using a broken power-law model. 



completed with radi o and IR data take n from NED, and the 
EGRET upper limit dFichtel et al. 1994 . 

We attempt to describe the two states in terms of changes 
in the physical parameters of the source. For both, we model 
the source as discussed in detail in Ghisellini et al. (2002). 
The model assumes synchrotron and Inverse Compton emis- 
sion from a spherical region in the jet, of radius R, at a distance 
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Fig. 2. Spectral Energy Distributions and models (solid lines) 
for Swift (red) and XMM-Newton+lSGRl (blue) data. The fig- 
ure shows in green arch ival radio and IR da ta (from NED) and 
the EGRET upper limit ( Fichtel et al. 1994 . The gray line rep- 
resents the GLAST sensitivity curve. The dashed line shows the 
input spectrum from the disk and the corona. 

Table 3. SED Model Parameters. See text for more details. 



Parameter 


Swift 


EPIC+ISGRI 


fl(10 15 cm) 


15 


15 


/i(10 15 cm) 


150 


150 


©(deg) 


3 


3 




15 


10 


L^lO^ergs- 1 ) 


8 


8 


y break 


5 


5 


/max 


1000 


1000 


B(G) 


3.3 


3.3 


W10 45 ergs-') 


150 


150 


R BLR (W l5 cm) 


1300 


1300 



h from the disk, with bulk Lorentz factor T and orientation 
with respect to the line of sight. The blob is populated by a 
tangled and homogeneous magnetic field B and by relativistic 
electrons, injected with total power L m j and a broken power- 
law energy distribution described by in [I, y break}, and y~ 3 - 2 
up to Jmax- The model accounts for radiative losses and deter- 
mines the cooling energy y coo \ of particles after a dynamical 
time R/c (see e.g., Maraschi et al. 2008). The primary emission 
is provided by a disk and a tenuous corona with a power-law 
spectrum of slope a = 1 and cut-off at 150 keV. The disk spec- 
trum is described by a black-body peaking at vo ~ 10 15 Hz and 
with power Ldisk, derived from the optical/UV spectral data. For 
the XMM-Newton+INTEGRAL state, we can place only weak 
constraints on the disk/corona emission, because of the lack of 
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data from the optical monitors of both satellites. However, the 
upper limit of XMM-Newton/OM in the UV band and the al- 
most constant archival optical data suggest that the soft photon 
input for the Inverse Compton regime does not vary signifi- 
cantly in the two states. The Broad Line Region (BLR), lo- 
cated at Rblr from the central black hole (see e.g., Ghisellini 
et al. 1996), reflects the pri mary emission toward the jet with 
luminosity Lblr - O.lLdisk dKaspi et al. 20051) . 

We attempted to reproduce both states by varying as small 
a number of parameters as possible. This choice was also mo- 
tivated by the limited energy range covered by the available 
simultaneous data. The absence of information above hundreds 
of keV makes the high energy electron distribution poorly de- 
termined. The synchrotron regime is also poorly constrained 
because of the lack of data. This is, however, unimportant for 
the SED modeling, since the radio-IR emission includes con- 
tributions from the outer regions of the jet and, therefore, is not 
considered in this model. 

Input model parameters for both states are given in Table[3] 
The overall SED of PKS 2 149 - 306 is typical of FSRQs at high 
redshift, with the Inverse Compton mechanism dominating the 
spectrum from X - to y-rays and in acco rdance with the standard 
blazar sequence dMaraschi et al. 2008). 

The two different SEDs can be explained in terms of a 
change only in the jet bulk Lorentz factor, being F = 15 for 
the "high-flux" state and F = 10 for the "low-flux" state. The 
different Lorentz factors of the emitting region can be inter- 
preted in terms of the inner-shock model scenario (e.g., Spada 
et al. 2001), which assumes that subsequent colliding shells at 
different speed produce a new radiating shell with intermedi- 
ate F. The variation in the F Lorentz factor accounts for the 
change in the total power of the source and for the spectral 
changes in the Inverse Compton regime, observed in the SEDs 
for the two states. In the "low-flux state" (found in XMM- 
Newton and INTEGRAL data), the "low" T drives a "cooling 
branch" in the electron distribution with y coo i -12 and a slope 
2.2 in [jpeak, 7a>oi] - This cooling branch does not develop in the 
"high-flux" state, since we found that y coo i ~ y break- 

The present data quality and the partial simultaneous spec- 
tral coverage limit a careful inspection of all full model input 
parameter range. However, it is remarkable that the two spectral 
states are described by changing only the jet bulk Lorentz fac- 
tor; this parameter accounts for both the flux variation and the 
spectral change of the Inverse Compton spectra. More complex 
solutions, changing more than one parameter, are not excluded, 
although we propose the simplest and physically consistent ex- 
planation of the state variation for the source. 



4. Conclusions 

We have investigated the long-term behavior of the blazar 
PKS 2149 - 306 (z = 2.345). All XMM-Newton, Swift, and 
INTEGRAL publicly available data have been reprocessed and 
analyzed with the latest software versions and calibration data 
bases. The SED, compiled with the data analyzed in the present 
work complemented with those available in the archives, is typ- 
ical of FSRQs at the brightest end of the blazar sequence. 



Swift observations suggest a hint of a roll-off below ~ 1 keV 
emerging on ~ 3 -day timescale, which can be explained as 
the natural spectral break due to the Inverse Compton onset. 
However, the present data quality does not allow us to distin- 
guish firmly between the absorption/curvature dichotomy, and 
further multiband and simultaneous observations are necessary. 

In the 60 - 300 keV energy band (in the blazar frame), the 
flux levels derived from ISGRI and BAT data differ by more 
than a factor of two, without spectral slope variation. From 
the broad-band spectral analysis, we can identify two differ- 
ent states: one found in data until 2004 {XMM-Newton and 
INTEGRAL), and another referring to a Swift observation per- 
formed in 2005. The SED modeling shows that the two states 
can be reproduced by changing only the bulk Lorentz factor of 
the jet. 
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